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Abstract 

The Bak-Sneppen model of self-organized biological evolution 
of an infinite ecosystem of randomly interacting species is repre- 
sented in terms of an infinite set of variables which can be consid- 
ered as an analog to the set of integrals of motion of completely 
integrable system. Each of this variables remains to be constant 
but its influence on the evolution process is restricted in time and 
after definite moment its value is excluded from description of the 
system dynamics. 

1 Introduction 

Investigations of complex dynamical systems play an important role in 
the modern theoretical and mathematical physics. In the recent years, 
the essential achievements have been arrived in studies of two dimen- 
sional field theoretical models. A large class of the exactly solvable and 
completely integrable models has been found [1]. It enabled one to un- 
derstand the laws of the soliton dynamics creating many non-ordinary 
phenomena in several physical systems. 

The non-trivial completely integrable models like sine-Gordon or non- 
linear Schrodinger equations are characterized by the infinite set of the 
conservation laws [1]. If the dynamics of such a system is considered 
in the framework of Hamilton's formalism and the integrals of motions 
are chosen as momenta, then the evolution in such a description appears 
to be similar to one of the mechanical system of non-interacting mate- 
rial points: the momenta are constant and the coordinates are linear 
functions of time. 

Many theoretical physicists hope that there are universal types of dy- 
namics of natural processes, which can be classified. A possible theoreti- 
cal approach for constructing such a classification could be the following. 
For the model under consideration one needs to find its "principal" rep- 
resentation, i.e. to formulate the dynamics of the system in terms of 
"principal" variables for which the dynamics is as simple as possible. If 



the principal representations of two models are the same, it is naturally 
to consider these models as representatives of the same type of dynamics. 
From this point of view the system of free particles is the representative 
of the class of all the completely integrable systems. 

If the construction of such a typology could be possible it would be 
very useful for understanding the most essential features in dynamics of 
complex natural systems. However, to find what we call the principal 
representation is not an easy problem for complex dynamical system. 
Probably, the set of ones for which this problem can be solved is restricted 
but it seems to be important to reveal the non-trivial systems of such a 
kind. 

In the present paper the proposed by Bak and Sneppen model (BSM) 
of biological evolution [2] is considered. It is a dynamical system describ- 
ing the ecosystem evolution as mutation and natural selection of inter- 
acting species. The most essential property of dynamics in this model 
is the self-organized criticality (SOC). The main characteristic of the 
SOC dynamics is that the system evolves to a critical state without fine 
tuning of its parameters. The SOC type in the BSM has the main spe- 
cific features of real biological evolution considered in the framework of 
the Gould-Eldridge "punctuated equilibrium" conception [3]. One can 
hope that this model represent an important universal type of critical 
dynamics of avalanche- like processes, realizing in many systems [4]. 

The simplest version of the BSM with random interaction structure 
[5] appeared to be exact solvable both in the thermodynamical limit and 
for finite system [6], [7], [8]. For infinite system it is possible to construct 
its principal representation. It is done in this paper. The dynamics in 
terms of principal variables looks as follows: each variable conserves its 
initial value the definite ( different for several variables ) length of time 
and after that becomes to be equal to zero. Thus the principal variables 
in this model are not exact integral of motion but in the process of the 
system evolution they transform (to zero) only once. For an alternative 
set of principal variables, one of them is lost at each time step and the 
others are renumbered. 

2 Formulation of the model 

In the framework of the BSM the biological evolution of ecosystem is 
described as follows [1]. The state of the ecosystem of N species is 
characterized by a set {x±, x^} of N number, < x$ < 1. The number 
Xi represents the barrier of the z-th species toward further evolution. 
Initially, each Xj is set to a randomly chosen value. At each time step 
the barrier x^ with minimal value and K — 1 other barriers are replaced 
by K new random numbers. In the random neighbor model (RNM) [5] 
which will be considered in this paper the K — 1 replaced non-minimal 
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barriers are chosen at random. In the local or nearest neighbor model 
(LM) these are the barriers of the nearest neighbors to the species with 
minimal barrier. For each species in the LM the nearest neighbors are 
assumed to be defined. 

For the LM the most of results are obtained by numerical experiments 
or in the framework of mean field approximation. The RNM is more 
convenient for analytical studies. The master equations obtained in [6] for 
RNM are very useful for this aim. These equations appeared to be exact 
solvable. The stationary solution was found in [6]. The time dependent 
solutions were obtained in [7] (infinite system), [8] (finite system). 

In this paper we construct the principal representation of the RNM 
master equation for infinite system. The obtained in [6] master equations 
for the RNM are of the form: 

P n {t + 1) = E ° l J( K - 0A'(1 - A)*-'P n _ {+ i(t) + 

1=0 

+8(K - n)QA n (l - X) K - n P (t) (1) 

Here, P n (t) is the probability that n is the number of barriers having 
values less than a fixed parameter A at the time t; < n < N, < A < 1, 
t > 0. For compact writing we used the threshold function 6(n) defined 
as follows 

9{n) = 1 for n > and 9(n) = for n < 0. 

The initial probability distribution P n (0) is supposed to be given. In 
virtue of the definition of P n (t), 

oo 

Pn{t)>0, £P„(*) = 1. (2) 
n=0 

Making summation in ([[]) over n it is easy to prove that 

N N 

£p„(f+i) = x;p B (f). (3) 

n=0 n=0 

In virtue of ([$]) if P n (0) are chosen in such a way that equalities (0) are 
fulfilled for t = 0, then for the solution of (0) it is the case for t > too. 
For analysis of (|l|) it is convenient to introduce the generating function 
q(z,u): 

oo 

q(z,u)= £ P„(f)*V. (4) 

t,n=0 

In virtue of (13) the function q(z,t) is analytical in z and u for \z\ < 1, 
\u\ < 1. The master equations (1) can be rewritten for the generating 
function q(z, u) as follows: 

(z-u(l-X + \z) K )q(z,u) = (z-l)u{l-X + Xz) K q(0,u) + zq{z,0). (5) 
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In d5|) the function 

N 



q(z,0) = J2Pn(0) 



z 



n=0 



is assumed to be given. If the generating function q(z, u) (P is known, 
P n (t) can be obtained as 



Qn+t 

dx n du t 



q(z, u)dzdu 

*=«=0 | 2 f =£1 \ u \ =t2 



z n+l u t+l 



where e\ < 1 and e 2 < 1. 



3 Principal representation 

Let us denote a(u) the analytical in u = solution of the algebraic 
equation 

a(u) - u{l + X(a{u) - 1)) K = 0. (6) 

It has the form 

a{u) = u{\ - \) K + Ku 2 X(l - X) K ^ + . . . . 

This series converges and \a(u)\ < 1 if \u\ < u and parameter u is 
chosen small enough. For K = 2 a(u) can be written as follows 



1 - 2A(1 - X)u - J(l - 4A(1 - X)u) 
a = a(u) = — I . 

Let us define the function j3(z): 

M - (1 . A Ur < 7 > 

It is analytical and \(3{z)\ < 1 if \z\ is small enough. In virtue of defini- 
tions (§, © 

a(/3(z)) = 2, /?(a(u)) = u. (8) 
Now, we define the function d(y,u) 

1 - a(y) 

which is analytical in y and u in the neighborhood of y = and it = 0: 

oo 

d(y,u)= £ C n (*)2/V. (10) 

n,t=0 



Substituting in (^) y — (3(z) and taking unto account (||) we obtain 

q(z,u) = (l-z)d((3(z),u). (11) 



The functions C n (t) denned by (|10D we consider as new dynamical 
variables of the RNM. The equalities (|9|), ( |TT| ) are the compact formulas 
of variable transformations from P n (t) to C n (t) and backward. They are 
not dependent in an evident way on the time what is seen if (^), ( |TTD is 
written in the more detiled form: 

oo oo oo oo n>(y) n 

E p n(t> n = (1 - *) E Cn(t)P(zr, E C n (t)y n = E P ^)t^7J7- V 

n=0 n=0 n=0 n=0 1 "li/J 

By substitution ([□]) in (Q) the following equation is obtained 

(/3(z) - u)d((3(z),u) = /3(z)d(j3(z),0) - ud(0, u). (12) 



Setting in fll2|) z = a(y) we have 

(2/ - u) = 0) - ud(0, it). (13) 

It follows from ( |13"D that 



■(y - u)d(y,u) 



= for n > 0, t > 0. 

y=u=0 



QyU+lQyt+1 

what looks in terms of the functions C n (t) as 

C n (t + 1) = C n+1 (t) for n > 0, t > 0. (14) 
The initial conditions 

C„(0) = c n (15) 

are dinned by q(z, 0): 

dM = ± Cny » = pm. (16) 

The equations (1^) has the simple solution: 

C n (t) = c n+t . (17) 

The variables C n (t) can be considered as the principal ones for the RNM. 
Principal representation of its dynamics is defined by equations (PI), 
with initial conditions ([15|), (|16|). Setting in (p~3|) it = y we see that 
0) = d(0, y) and 

./ \ yd(y,0)-ud(u,Q) 

d{y,u) = . (18) 

y-u 
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Setting y = /3(z) in ( pT8[ ) and taking into account (|9]), (|TTD we obtain the 
solution of the master equation (||). 

_ z(l - a(tt))g(z, 0) + u(z -!)(! + A(g - l)) K q(a(u), u) 
q ^ U >~ (l-a(u))[z-u(l + \(z-l)) K \ 

If time dependent transformations of variables are allowed, an alter- 
native principal representation can be constructed. Let us define the 
variables S n (t) as follows 

S n (t) = C n - t (t) for t < n, and S n (t) = for t > n; n > 0, t > 0. 
The backward variable transformation can be written in the form 



Cjt) = Z n+t (t), n>0, t>0. 
It follows from ( |17D that 

S n (t) = c n 6(n - t), 

i.e. the variable S n (t) conserves its initial value S n (0) = c n until t < n 
and S n (t) = from the moment t — n + 1. 



4 Conclusion 

We have obtained the following results. In terms of variables chosen in a 
special way the SOC dynamics in the BSM can be described by very sim- 
ple equations (0), fll4p- These variables are expressed straightforwardly 



in terms of infinite set of the constants defined explicitly by the initial 
conditions. On each time step one of these constants is forgotten, i.e its 
value does not influence the further stages of the system evolution. The 
consequent loose of the information about initial state is all what hap- 
pens in the self-organization process described by RNM. The system of 
such a kind could be called completely integrable dissipative system. It 
would be interesting to understand how robust this dynamics is. Is it the 
inherent property of the SOC processes or an artifact of the considered 
model? 
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